home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / dpbcon.z / dpbcon
Text File  |  1996-03-14  |  3KB  |  133 lines
  1.  
  2.  
  3.  
  4. DDDDPPPPBBBBCCCCOOOONNNN((((3333FFFF))))                                                          DDDDPPPPBBBBCCCCOOOONNNN((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      DPBCON - estimate the reciprocal of the condition number (in the 1-norm)
  10.      of a real symmetric positive definite band matrix using the Cholesky
  11.      factorization A = U**T*U or A = L*L**T computed by DPBTRF
  12.  
  13. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  14.      SUBROUTINE DPBCON( UPLO, N, KD, AB, LDAB, ANORM, RCOND, WORK, IWORK, INFO
  15.                         )
  16.  
  17.          CHARACTER      UPLO
  18.  
  19.          INTEGER        INFO, KD, LDAB, N
  20.  
  21.          DOUBLE         PRECISION ANORM, RCOND
  22.  
  23.          INTEGER        IWORK( * )
  24.  
  25.          DOUBLE         PRECISION AB( LDAB, * ), WORK( * )
  26.  
  27. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  28.      DPBCON estimates the reciprocal of the condition number (in the 1-norm)
  29.      of a real symmetric positive definite band matrix using the Cholesky
  30.      factorization A = U**T*U or A = L*L**T computed by DPBTRF.
  31.  
  32.      An estimate is obtained for norm(inv(A)), and the reciprocal of the
  33.      condition number is computed as RCOND = 1 / (ANORM * norm(inv(A))).
  34.  
  35.  
  36. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  37.      UPLO    (input) CHARACTER*1
  38.              = 'U':  Upper triangular factor stored in AB;
  39.              = 'L':  Lower triangular factor stored in AB.
  40.  
  41.      N       (input) INTEGER
  42.              The order of the matrix A.  N >= 0.
  43.  
  44.      KD      (input) INTEGER
  45.              The number of superdiagonals of the matrix A if UPLO = 'U', or
  46.              the number of subdiagonals if UPLO = 'L'.  KD >= 0.
  47.  
  48.      AB      (input) DOUBLE PRECISION array, dimension (LDAB,N)
  49.              The triangular factor U or L from the Cholesky factorization A =
  50.              U**T*U or A = L*L**T of the band matrix A, stored in the first
  51.              KD+1 rows of the array.  The j-th column of U or L is stored in
  52.              the j-th column of the array AB as follows:  if UPLO ='U',
  53.              AB(kd+1+i-j,j) = U(i,j) for max(1,j-kd)<=i<=j; if UPLO ='L',
  54.              AB(1+i-j,j)    = L(i,j) for j<=i<=min(n,j+kd).
  55.  
  56.      LDAB    (input) INTEGER
  57.              The leading dimension of the array AB.  LDAB >= KD+1.
  58.  
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. DDDDPPPPBBBBCCCCOOOONNNN((((3333FFFF))))                                                          DDDDPPPPBBBBCCCCOOOONNNN((((3333FFFF))))
  71.  
  72.  
  73.  
  74.      ANORM   (input) DOUBLE PRECISION
  75.              The 1-norm (or infinity-norm) of the symmetric band matrix A.
  76.  
  77.      RCOND   (output) DOUBLE PRECISION
  78.              The reciprocal of the condition number of the matrix A, computed
  79.              as RCOND = 1/(ANORM * AINVNM), where AINVNM is an estimate of the
  80.              1-norm of inv(A) computed in this routine.
  81.  
  82.      WORK    (workspace) DOUBLE PRECISION array, dimension (3*N)
  83.  
  84.      IWORK   (workspace) INTEGER array, dimension (N)
  85.  
  86.      INFO    (output) INTEGER
  87.              = 0:  successful exit
  88.              < 0:  if INFO = -i, the i-th argument had an illegal value
  89.  
  90.  
  91.  
  92.  
  93.  
  94.  
  95.  
  96.  
  97.  
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.